Chaotic dephasing in a double-slit scattering experiment 
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We design a computational experiment in which a quantum particle tunnels into a billiard of variable 
shape and scatters out of it through a double-slit opening on the billiard's base. The interference 
patterns produced by the scattered probability currents for a range of energies are investigated in 
relation to the billiard's geometry which is connected to its classical integrability. Four billiards 
with hierarchical integrability levels are considered: integrable, pseudo-integrable, weak-mixing and 
strongly chaotic. In agreement with the earlier result by Casati and Prosen (Phys. Rev. A 72, 
032111 (2005)), we find the billiard's integrability to have a crucial influence on the properties of the 
interference patterns. In the integrable case most experiment outcomes are found to be consistent 
with the constructive interference occurring in the usual double-slit experiment. In contrast to 
this, non-integrable billiards typically display asymmetric interference patterns of smaller visibility 
characterized by weakly correlated wave function values at the two slits. Our findings indicate an 
intrinsic connection between the classical integrability and the quantum dephasing, responsible for 
the destruction of interference. 



The chaotic phenomena in dynamical systems 
have been extensively studied and w^ell under- 
stood over the last fev^ decades. Yet, the quan- 
tum manifestations of classical chaos, despite 
largely investigated, are still not entirely clear. 
A possible bridge between the two is provided 
by the mechanism of quantum dynamical de- 
phasing, occurring when a coherent quantum 
state undergoes a classically chaotic dynamics. 
In this work we report the results of a compu- 
tational experiment that models scattering of a 
quantum particle through billiards whose classi- 
cal dynamical properties depend on their shape. 
By relating the scattering outcomes with the 
integrability of the billiards (scatterers), we es- 
tablish a clear connection between the classical 
chaos and the appearance of quantum dephas- 
ing. In particular, we find that a large disper- 
sion of interference visibilities (contrasts) with 
respect to the energy sampling suggests a signa- 
ture of classical chaos. Our results are in agree- 
ment with the previous computational and real 
experiments, and can themselves be obtained 
in a real laboratory experiment. 



1 Introduction 

The precise relationship between quantum and classi- 
cal dynamics to this date still remains somewhat elu- 
sive, at both theoretical and experimental level. In 
particular, the understanding of quantum counterpart 
of macroscopic chaotic behavior is far from complete 
[21 [3] . As pointed out in [1] , there are two crucial prop- 
erties of quantum motion that are in a direct contrast 
with common characteristics of classical chaos, namely: 
(i) discrete energy-spectrum of bound quantum sys- 
tems, rendering quantum dynamics quasi-periodic and 
hence decomposable in a finite Fourier-spectrum; (m) 
inherent stability of quantum dynamics, as a result of 
which initial errors always propagate only linearly with 
time. 

Nevertheless, various general manifestations of clas- 
sically chaotic phenomena in quantum dynamics have 
been widely studied over the last few decades [21 [3l 
m [5]. Recently, some of the long-standing theoreti- 
cal results in this direction received their experimental 
confirmation, such as the chaotic nature of quantum 
motion of the kicked top [6]. 

On the other hand, the phenomena of quantum 
decoherence provides a model to understand how the 
macroscopic classical dynamics arises from the micro- 
scopic quantum dynamics [?]• Dynamical decoherence, 
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specific form of quantum decoherence, was pro- 
posed in the context of studying the dynamical behav- 
ior of a quantum analog of a classically chaotic sys- 
tem [H O [9]. In particular, it has been recognized 
that the classically chaotic dynamics is related to the 
mechanism of quantum dephasing, responsible for the 
destruction of the phase-coherence of the initial wave 
function. 

Time-dependent experiment. In their recent work 
[T], Casati and Prosen designed a computational exper- 
iment in order to expose the mechanism of quantum 
dephasing arising due to (classically) chaotic dynam- 
ics. The time-dependent Schrodinger equation: 



was solved numerically for a quantum particle moving 
freely inside the two-dimensional domain indicated by 
full line in Fig.[TJ The particle was initially prepared in 
a minimum uncertainty coherent state and placed in- 
side the upper triangular region (billiard domain) with 
a given initial momentum. Two small slits were con- 
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Figure 1: (from [T]) Time-dependent experiment: bil- 
liard with a variable geometry connected to the radiat- 
ing domain through a double-slit opening. An initially 
coherent state freely evolves within the domain, creat- 
ing interference patterns on the screen. 

necting the billiard domain with the lower radiating 
domain, allowing the particle to leak out in form of 
two interfering probability currents during the time- 
evolution. The particle's wave number k defining the 
de Broglie wavelength A = 2Tr/k was chosen in a way 
to have A much bigger than the width of the slits, and 
much smaller than the distance between them. The 
interference patterns, created by time- integrating the 
perpendicular component of the probability current on 



the horizontal screen below, were observed. The exper- 
iment went on until the probability to find the particle 
within the billiard domain was vanishingly small. 

The experiment was repeated for two different sha- 
pes of billiard domains: the right-angle isosceles trian- 
gle, and the Sinai-type (defocusing) billiard obtained 
by replacing this triangle's hypotenuse with a circu- 
lar arc (dashed line in Fig.[T]). These two billiards 
have drastically different classical dynamical proper- 
ties: while the former generates regular integrable mo- 
tion, the latter is known to be a strongly chaotic K- 
system (strong mixing) [10]. Fundamentally different 
results were found for the integrated interference pat- 
terns, as shown in Fig.[2j In the integrable case the 
usual double-slit fringes characterizing constructive in- 
terference were recovered, whereas in the chaotic case 
the fringes were entirely absent and only a bell-shaped 
unimodal intensity distribution was found. In particu- 
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Figure 2: (from [1], color online) Results of the time- 
dependent experiment: while the integrated interfer- 
ence pattern for the case of integrable billiard shows 
constructive double-slit fringes, the pattern displayed 
in the chaotic case appears to be identical to the su- 
perposition of two uncorrelated sources. 

lar, the chaotic billiard's profile was essentially iden- 
tical to the situation of two slits emitting indepen- 
dent probability currents. While the classically reg- 
ular /integrable dynamics preserved the initial phase- 
coherence of the particle, the chaotic dynamics pro- 
duced pseudo-random phase shifting, ultimately gener- 
ating "uncorrelated" output currents. Classically chaotic 
time-evolution induced the quantum dephasing pro- 
cess, responsible for the destruction of particle's phase- 
coherence, and thus causing the disappearance of the 
expected double-slit fringes. 

This result has been recently confirmed in a real ex- 
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periment investigating wave chaos [TT]. Furthermore, 
similar computational experiments have confirmed an 
intrinsic relationship between classical chaos and quan- 
tum dephasing in a wide range of examples |8|, while 
the important role of the geometric symmetry has been 
pointed out in [9]. 

In this work we present the results of a new compu- 
tational experiment, constructed along the similar lines 
as the time-dependent experiment described above, but 
now considering the energy domain instead of the time 
domain. We consider stationary scattering of quantum 
particles with different energies through the billiards of 
variable geometries, designed to include a double-slit 
opening. The interference patterns are observed on the 
screen over a wide energy-interval, and statistically in- 
vestigated in relation to the billiards' geometries. As 
we show, the findings from time-dependent experiment 
fully extend to the scattering case, shedding new light 
on the mechanism of dynamical dephasing. Moreover, 
two new billiards with the particular non-integrable 
geometries, whose dynamics interpolates between in- 
tegrable and strongly chaotic, are added in order to 
study the appearance of dephasing in relation to the 
ergodic hierarchy of dynamical systems. In opposition 
to [9], we avoid the presence of any kind of geometric 
symmetries of our billiards, and focus exclusively on 
their integrability properties. 

For each energy value we obtain an interference pat- 
tern whose properties are essentially given by the dif- 
ference and ratio of the wave function's phases and 
magnitudes, respectively, at the positions of the two 
slits. We specifically consider the distribution of the 
interference visibilities, given as the range of ratios of 
the wave function intensities. Our main finding is that 
in the case of chaotic classical dynamics, this distri- 
bution is much wider than in the case of integrable 
dynamics: a "typical" interference pattern produced 
by each of the non-integrable billiards is much more 
likely have lower visibility than the one created by the 
integrable billiard. 

This paper is organized as follows: in the next Sec- 
tion we describe in detail the design of our computa- 
tional experiment and the parameters used. In Section 
[3] we expose the results found, and finally provide a 
discussion and conclusions in Section [H 

2 Scattering experiment 

In this Section we explain the construction of our scat- 
tering apparatus, provide the parameter values used in 
the experiment, and outline the computational details. 



We consider a isosceles right triangle billiard with 
the double-slit opening on its base (as the one used in 
time-dependent experiment), shown in Fig. [3^. We re- 
place its left side with a potential tunneling barrier of 
height Vbarrier shown in red, and use it to connect the 
billiard to the input channel having the width equal to 
that triangle's side. All other billiard borders shown 
in black are assumed to be infinite potential barriers. 
The billiard 's base containing the slits is connected to 
the output channel having its width equal to it. In- 
side the output channel at a distance dgcreen from the 
base, we put the screen where the interference patterns 
are observed and recorded. Thus constructed scatter- 
ing system consists of three parts: input channel used 
to inject the particle into the triangular billiard (illus- 
trated as cyan arrows), the billiard itself used as the 
scatterer, and the output channel where the interfer- 
ence patterns are observed. The scattering experiment 
is done for four different billiards, out of which we have 
just described the first. Other billiards are designed 
by modifying this basic triangular billiard, and named 
through their classical dynamical properties. In par- 
ticular: 

• integrable - (just described) isosceles right-angle 
triangle with angles f,f,f, shown in Fig. [3^. 
Classically, it generates exclusively regular, pe- 
riodic or quasi-periodic motion. 

• pseudo-integrable - right-angle triangle with an- 
gles 5 1 f 5 classically a pseudo-integrable sys- 
tem |12] . shown in Fig. [3b. 

• weak-mixing - right-angle triangle with two an- 
gles irrational with vr, classically conjectured to 
be a weak-mixing system [13], shown in Fig.[3h. 
For angles we choose 7r ^~^ ,tt ^^^ , |. 

• chaotic - the billiard obtained by replacing the 
triangle's hypotenouse with a circular arc of ra- 
dius R, as done in the time-dependent experi- 
ment, known to be a strongly chaotic K-system 
[10], shown in Fig.[3li. 

A quantum particle in the form of a plane wave with 
definite energy E = is injected into the billiard 
through the input channel by tunneling the potential 
barrier, and is then scattered out through the slits. 
We model this process by numerically solving the "on- 
shell" time- independent scattering problem for the two- 
dimensional Schrodinger equation 
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Figure 3: (color online) Scattering experiment set-up involving four billiards with different classical dynamical 
properties: (a) integrable, (b) pseudo-integrable, (c) weak-mixing, and (d) chaotic. Billiards are attached 
to an input channel of unitary width and separated from it by a potential tunneling barrier that replaces 
that billiard's side. Base of each billiard contains a double-slit opening, leading to the output channel where 
the interference patterns are recorded on the screen at a distance dgcreen- Particle's wavelengths are shown 
for a comparison with slits width s, distance between the slits dg, and width of the circular segment v (for 
chaotic billiard in (d)). Arrows illustrate the momentum of the incoming particles (plane waves) that impact 
orthogonally on the potential barrier. We take the plane wave direction to be x-coordinate, and the direction 
of the scattered currents to be y-coordinate. In all the figures to follow, unless noted otherwise, we use this 
scheme for representing the results for four billiards. 



with scattering boundary conditions. The scattering 
wave function 'I' consists of three spatially disjoint parts: 

• 'I'l in the input channel - incoming plane wave 

• ^2 inside the billiard - scattering state 

• ^3 in the output channel - scattered waves 

The final solution ^ = ^i©^'2®^3 is obtained by im- 
posing the spatial continuity condition between and 
^2 5 and between 'I'2 and ^3, along with the boundary 
conditions given by the billiard's geometry (infinite po- 
tential barriers) and the incoming particle. While us- 
ing as input, we will investigate the properties of 
^2) and in particular focus on ^'3 which generates the 
interference patterns related to the billiard's integra- 
bility properties and the particle's energy E. 

Majority of the probability current is immediately 
reflected back from the potential barrier and does not 



undergo the scattering process. However, if penetrates 
the tunneling barrier, the particle is much more likely 
to exit the billiard through the slits than back through 
the input channel (penetrating the barrier again) . The 
presence of the barrier drastically reduces the intensity 
of the probability currents exiting the billiard and in- 
terfering, but on the other hand allows the penetrating 
particle to "feel" the billiard's geometry much better. 
The barrier's height is kept equal for all billiards. In 
opposition to the time-dependent experiment, we are 
now not integrating the recorded interference patterns 
(in energy), but instead examine their properties sta- 
tistically over the entire energy-interval. Also, we do 
not use any absorbing potential at the edges of the 
output channel. We took the input channel to have 
the maximum possible width (that of the billiard), in 
order to obtain the best focus in particle's direction 
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(momentum) while allowing for maximal indetermi- 
nacy in position. This way we optimize the directional- 
ity of the plane wave's impact on the barrier, assuring 
that the scattering outcome as best possible reports 
about the properties of the billiard. Besides the bar- 
rier, the strength of probability currents also depends 
on the eigen-energies of the closed billiard, that act 
as resonances in the scattering process. The transport 
through the scattering apparatus is enhanced when the 
particle's energy corresponds to one of the scatterer's 
(billiard's) resonant energies. However, we will con- 
sider an ensemble of all the recorded interference pat- 
terns with equal statistical weights, regardless of the 
particles' energies. 

By considering the four billiards we investigate the 
scattering process in relation to a hierarchy of systems 
having different classical "chaoticity" levels. Thus, we 
obtain a more detailed relationship between the classi- 
cal dynamical properties and the appearance/strength 
of the quantum dephasing. For all billiards we set the 
width of the input channel (and hence the length of the 
billiard's potential barrier side) to be 1, therefore fixing 
the billiards' areas. This implies that the billiards will 
have slightly different areas, while pseudo-integrable 
and weak-mixing billiards will also have their base side 
somewhat shorter (and hence their output channels 
somewhat narrower). The chaotic billiard is designed 
with a fixed circular arc replacing the hypotenouse, 
and with the base equal to 1 (that is to say, by cut- 
ting a circular segment away from the isosceles triangle 
billiard). 

The parameter values used. We consider the sys- 
tem of units such that not only the width of the bil- 
liard, but also the Planck's constant and the particle's 
mass equal unity, h = m = 1. We seek to have slits 
width s much smaller than the particle's wavelength A, 
which itself has to be much smaller than the distance 
between the slits d^: 

s <C A <C 4 (3) 

Furthermore, for the chaotic billiard we take the width 
of the cut circular segment to be larger (but compa- 
rable) to the particle's wavelength A (cf. Fig.[3]i): 



in order for "chaoticity" of this billiard to be "visible" 
by the quantum particle. Accordingly, we set the pa- 
rameters and dimensions used in our experiments as 
follows: 

* the lengths of the billiards' sides - integrable: 1, 1, 1.41 



pseudo-integrable: 1, 0.727, 1.236; weak-mixing: 1,0.684, 
1.212; chaotic: 1,1,1.428. 

* the billiard's areas A - integrable: 0.5; pseudo-integrable: 
0.363; weak-mixing: 0.342; chaotic: 0.414. 

* curvature radius for chaotic billiard: R = 2.786, the 
width of the cut circular segment: ly = 0.091. 

* width of the slits s - integrable and chaotic: 0.017, 
pseudo-integrable and weak-mixing: 0.016. 

* distance between the slits dg ~ integrable and chaotic: 
0.287, pseudo-integrable and weak-mixing: 0.257. 

* the energy-interval of incoming particles was mea- 
sured in terms of particle's wave number k = V^E, 
and set from kmin = 85 {Xmax = 0.074, Emin = 3612), 
to kmax = 94 {X^in = 0M7,E^ax = 4418). This 
choice optimizes the condition Eq. Q, along with the 
condition A < stated above. 

* the fc-step between two particles of consecutive k- 
values: 6k = 0.001. We made 9000 experiments for 
each billiard, with particles uniformly spaced in k, cov- 
ering the interval of k mentioned above. 

* the height of the potential tunneling barrier: Vhamer = 

7200, corresponding to kbarrier = V^Warrier = 120. 

This value was set by maximizing the dephasing effects 
while still having measurable probability currents (see 
Section |4]). The barrier's energy is thus roughly twice 
bigger than average particles' energies. 

* distance from the billiard's base at which the inter- 
ference patterns are observed: dscreen — 

0.12. 

Therefore, in our system of units the Eq. ([2]) reduces 
to: 

We numerically solved this equation (with appropriate 
boundary conditions) for all fc-values and parameters 
as described above, and for all four billiards. 

The numerical implementation. The Eq. is 

solved numerically by discretizing the space where the 
wave function ^ is defined. The space within the bil- 

BiB I 1) 

liards is covered by a grid of — - mesh points, the 
potential barrier, input and output channels are each 
covered by a sequence of B points, while each of the 
slits is covered by S points. Each point is labeled in 
two dimensions as ipn,m, and represents the wave func- 
tion complex value at its coordinates. The Laplacian 
is discretized as follows: for a given point Vn,m (e.g. 
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inside the bihiard), the Eq. ^ reads: 

• {lpn+l,m + V'n,-l,m " 2'0n,m) + 

) + (5) 

where A^^ and Ay are the space discretization steps 
in horizontal and vertical direction. This way we con- 
struct a linear system of dimensionality = ^^^^^^ -)_ 
3B + 25 in the form 

A^l^ = ffo (6) 

where the A^-dimensional vector ip is the solution of 
the scattering problem containing the points '4>n,m, vec- 
tor ffo contains the inputs defining the incoming plane 
wave (and zeros for other points as in example of Eq. ([H}). 
while the N x N matrix A contains discretized two- 
dimensional Laplacian and the conversion coefficients 
for other parts of the scattering apparatus. The vector 
'ip includes B values defining the wave function inside 
the output channel, used to compute the interference 
pattern. 

The grid points for the integrable billiard form a 
regular square lattice A^^ = A^. For other billiards 
we construct the appropriate grids in accordance with 
their geometries, so that all the boundary grid points 
lie exactly on the geometrical billiard's boundary. For 
the chaotic billiard we use the same 'trick' as in the 
time-dependent experiment: the spacing between the 
grid's horizontal rows Ay is made non-uniform in or- 
der to produce a smooth curvature modeling the cir- 
cular arc that replaces the hypotenouse. For pseudo- 
integrable and weak-mixing billiards, the spacing be- 
tween the grid's vertical columns A^ is reduced by a 
factor of tan ^ and tan7r ^~^ respectively, in order to 
model their geometries. In our experiments, we used 
the values oi B = 180 for all billiards (whereas we 
have carefully checked that numerics is robust against 
variation of B in this range), along with 5 = 3 for in- 
tegrable and chaotic, and 5 = 4 for pseudo-integrable 
and weak-mixing billiards (this explains slightly differ- 
ent values of slits distance and width for these billiards, 
as described in the previous paragraph). 

The total numerical error in the experiments was 
estimated by computing the sum of the probability flux 
through the slits (pout and the probability flux back- 
scattered through the input channel (barrier) (pbacki 
and comparing this sum with the probability flux of 
the incoming wave - normalized to 1. The total uni- 
tarity error given by \{(pout + (pback) -M was typically of 
the order of 10~^, and only in the vicinity of billiards' 
resonances went up to 10~^. 



3 Results 

In this Section we expose the results of the compu- 
tational experiments designed and carried out as de- 
scribed in the previous Section. We systematically 
compare the findings among four billiards. 

The scattering states and resonance profiles. We 

examine the billiard's scattering states given by the 
wave function ^'2 defined within the billiards. In Fig. [4] 
we show the color-plot of the densities |^'2p for all 
billiards for the wave number A: = 90. The wave func- 
tion density for the integrable billiard is clearly very 
regular, while the chaotic billiard's density appears 
to be rather irregular. Pseudo-integrable billiard dis- 
plays a density structurally similar to the integrable 
case, while the weak-mixing density seems closer to 
the chaotic one. The regularity of the wave function 
densities neatly follows the hierarchy of the billiards' 
classical integrability levels. The scattering wave func- 
tion's properties refiect the billiards' classical dynam- 
ical properties. The characteristic lengths of densities 
from Fig. [4] correspond to the incoming particle's wave- 
length. They do not vary substantially among the bil- 
liards, as A; = 27r/A is the same for all billiards shown 
in Fig.m The densities for different wave numbers 
k qualitatively preserve the regularity/irregularity of 
their structural patterns. 

We define the probability flux leaking through the 
slits (p 

out as the total probability current passing through 
any perpendicular cross-section Cout of the output chan- 
nel (such as the screen for instance): 

(pout = -Im / ^-^'^s 

which sensitively depends on the particle's energy, and 
its relationship with the billiard's resonant eigen-ener- 
gies. We investigate the values of the flux (pout as a 
function of the particle's wave number k, and report 
our findings in the Fig.[5j For all the four billiards the 
profile of (/)out-values displays peaks indicating the posi- 
tions of billiards' scattering resonances (roughly corre- 
sponding to the eigen-energies of the closed billiards). 
The density of peaks slightly varies among the billiards 
as a consequence of their different areas. This density 
is proportional to the billiard's area A, which is 
known as the Weyl's law |14j : 

dNr _ A ^ 
"dk ~ 2tt 

(for this reason we used particles equally spaced in 
wave number rather than in energy). Since the wave 
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Figure 4: (color online) Densities of the scattering states 1^*2 P inside the billiards for the A;- value of k = 90. 
The scale of |^'2p is shown on the right of each billiard. Big dots indicate the positions of the slits. Plots 
correspond to billiards as depicted in Fig.[3l 




Figure 5: Profiles of (/)o«t-value as s function of incoming particle's wave number k for all the four billiards, for 
the wave number interval from = 90 to = 92. Log-scale is used in (j)out display more clearly the structure 
of the resonances. Profiles correspond to billiards as depicted in Fig.[3l 
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number interval is the same for all billiards in Fig.[5l 
the ratio of number of peaks roughly follows the ratio 
of billiards' areas, i.e. 0.5 : 0.363 : 0.342 : 0.414. From 
this formula one can also estimate the total number of 
peaks in the given fc-interval. Note the different shapes 
of peaks: in the integrable case, immediately after 
the resonant fc-value, (pout sharply decays, whereas for 
other billiards the resonant peaks are relatively sym- 
metric. Also note that the densities of Fig. H] evaluated 
for k = 90 are on different scales (weak-mixing case 
in particular). This is a consequence of the particle 
with k = 90 producing different scattering outputs for 
different billiards: as clear from Fig.[5l for k = 90 the 
(/)out-value is biggest in the weak-mixing case. 

The interference patterns. In the reminder of this 
Section, we analyze in detail the interference properties 
of the scattered probability currents, comparing them 
among the billiards. The intensity I{x) of the vertical 
component of the probability current within the out- 
put channel, computed at the distance dgcreen from the 
billiard's base is: 



through both slits with equal intensities, which is di- 



I{x) 



Im 



(7) 



f'screen 



where x is the horizontal coordinate parameterizing the 
screen. For an interference pattern we intend the 
normalized intensity I{x): 



S(x) 



I{x) 



f I(x)dx 

J screen V ' 



This way we eliminate the pattern's dependence of the 
total output flux (poutik) and focus only on its interfer- 
ence properties, such as symmetry and visibility (con- 
trast ratio of intensity between adjacent minima and 
maxima) . Because of being a result of the interference 
between two probability currents, the patterns' proper- 
ties are linked to the wave function's values at the cen- 
ters of the (almost point-like) slits "^si and \I's2 (con- 
necting and ^s). Namely, we can approximately 
imagine the field ^3 to be a result of a simple inter- 
ference between two point sources, ^si and ^'<j2. In 
order to examine the patterns S(x) systematically, we 
introduce the quantities Ap and A^ as follows: 



A. 



2 Re In 



Im I In 



£l 



In- 



arg ^'si - arg 



(8) 



Ap reports the ratio of the wave function densities at 
the slits. Thus, Ap = means that the particle leaks 



rectly related to the interference visibility. A, 







implies the maximal visibility and hence the infinite 
contrast ratio in the far field (i.e. perfectly construc- 
tive interference corresponding to the usual text-book 
double-slit experiment). For increasing |Ap| visibility 
gradually worsens, and for |Ap| 1 the interference 
pattern is practically washed away. A^ gives the value 
of the phase-shift between the wave function values at 
the slits and directly controls the (a)symmetry of the 
interference pattern. = indicates a perfectly sym- 
metric S. Numerically, we compute the ^si,s2 by aver- 
aging the values for all the space-discretization points 
inside the slits (slits are very small anyway due to 
s <C A). 

In Fig. [6] we show the profiles of Ap and A^ as func- 
tions of the wave number k. Note a very different be- 
haviors for different billiards. The integrable case ex- 
hibits a smooth behavior with mostly Ap = A,^ = 0, in 
addition to sporadic jumps overlapping with the reso- 
nances. In contrast to this, the non-integrable billiards 
display a much more dispersed fluctuating behaviors 
that intensify with the billiards' "chaoticity" level. In 
particular, the Ap profile (although relatively smooth) 
shows large variations, mostly staying far from Ap = 
even for A;- values away from the resonances. In the in- 
tegrable case, the resonances are very narrow and the 
transport through the cavity is very small for most val- 
ues of k. This in particular corresponds to the regions 
where Ap and A^ are very small. Nevertheless, the 
variations of the interference visibility Ap are smaller 
in the integrable case than in the non-integrable cases 
even in the resonance regions (see also Fig. [9]). 

The pair of values (Ap,A(^) captures the essential 
properties of an interference pattern, quantifying its 
discrepancy from the constructive S. This is illustrated 
in Fig. [7] where we show the interference patterns for 
five characteristic examples of (Ap, A0)-values. The in- 
terference pattern shown in Fig.[71: having (Ap, A,^) = 
(0, 0) is produced by the integrable billiard, and it is 
essentially identical to the pattern obtained from the 
ideal double-slit experiment with perpendicularly im- 
pacting plane waves (with the appropriate particle's 
energy). Any other pattern with (Ap,A(^) 7^ (0,0) 
shows (at least) some departure from this ideal case. 
This is illustrated in Fig.[7] (a),(b),(d) and (e), where 
we show four typical examples of E-s generated by the 
chaotic billiard. For instance, the pattern in Fig.[7ti 
with Ap > is characterized by almost entire prob- 
ability fiux passing through only one slit. The pat- 
tern from Fig. [7^ is formed by the current exiting the 
billiard through one slit, in addition to a phase-shift 
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Figure 6: (color online) Profiles of A^-values in red and A^- values in blue (cf. Eq. ([8])) as functions of incoming 
particle's wave number for all the four billiards, for the A;-interval from k = 90 to k = 92. Profiles correspond 
to billiards as in Fig.[3l 





Figure 7: Five examples of different interference patterns in relation to their (Ap, A(^)-values. The case (c) 
is an example of a constructive interference pattern with Ap = IS.^ = 0, produced by the integrable billiard 
for k = 89.707. Other cases are examples of non-constructive and/or asymmetric patterns, as either Ap ^ 
and/or A^ / 0. They are generated by the chaotic bilhard, for: (a) k = 90.471, (b) k = 92.623, (d) k = 87.208 
and (e) k = 91.233. 



towards the other slit (Ap > 0, A^ < 0). This is 
also clear from Fig. [8] where we show two-dimensional 
interference of probability currents in the output chan- 



nel corresponding to the five situations from Fig. [71 It 
is only in the case from Fig. [8b that we have a fully 
symmetric two-dimensional interference characterized 
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Figure 8: (color online) Two-dimensional visualization of interference of two probability currents inside the 
output channel. The scale indicates the intensity of the vertical component of the current defined in Eq. ([7]). 
Five situations correspond to the five cases of interference patterns from Fig. [71 Thick line indicates the position 
of the screen. 



by both slits emitting equally strongly and in phase 
coherence. 

The interference properties of a specific pattern 
are basically independent of the billiard, and de- 
pend exclusively on the pattern's {Ap, A^)-valne. We 
verified this claim for vast range of patterns produced 
by different billiards with the same (Ap, A(^)-values. 
The pattern's properties do depend on the incoming 
particle's A;- value in the sense of total number of fringes. 
However, since we focus only on the contrast of T,{x), 
we disregard this dependence and consider all the pat- 
terns with the same (Ap, A^)-values to be identical 
(from the contrast /symmetry point of view). Although 
the patterns from Figs. [7]&[8] were obtained from inte- 
grable and chaotic billiards for different A;-values, they 
contain the properties of any billiard's S with those 
(Ap, A(^)-values. In the rest of this work we will iden- 
tify S-s only through their (Ap, A0)-values, and in 
particular, distinguish only between the constructive 
patterns defined by {Ap,A^) = (0,0) and other non- 
constructive patterns. 

The distribution of interference patterns. The 

remaining question we answer in the rest of this Sec- 
tion refers to the typical proportion of constructive 
versus non-constructive patterns, in relation to four 
billiards. After establishing a connection between a 
pattern's constructiveness properties and its (Ap, A^)- 



value, we statistically study the distributions of Ap 
and A(^. In Fig. [9] we show two-dimensional scatter 
plot distribution of Ap-A,^, done over the entire ex- 
amined /c-interval. In the integrable case the patterns 
are visibly concentrated around a single central peak at 
(Ap, A<^) = (0, 0). In opposition to this, non-integrable 
billiards display rather diffused and more uniform dis- 
tributions, without a clear central peak. While in the 
integrable case the constructive patterns are overwhelm- 
ingly present over a large /c-interval, other billiards 
exhibit more significant fractions of non-constructive 
patterns with large |Ap|, but also larger fraction of 
asymmetric patterns with large |A(^|. We illustrate 
this further by showing in Fig. [10] projected (marginal) 
distributions of Ap and A,^, simultaneously compar- 
ing all the four billiards. The distributions of Ap- 
values, P{Ap), shown on the left reports a clear dis- 
crepancy between the integrable billiard and the non- 
integrable billiards. P{Ap) in the integrable case ex- 
hibits a clear peak which exponentially decays away 
from Ap = 0. On the other hand, P(Ap) for other bil- 
liards show a gaussian-like profile with much larger dis- 
persion. The distributions of phase-shifts P(A(^) dis- 
play less prominent discrepancy among the billiards, 
although the same trend is visible. While statistics 
for both Ap and A^ show a distinction between inte- 
grable billiard and non-integrable billiards, the differ- 
ence among the non-integrable billiards themselves is 
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Figure 9: (color online) Two-dimensional distributions of (Ap, A0)-values obtained for all experiments done in 
the fc-interval from k = 85 to k = 94. The fraction of the interference patterns having a certain (Ap, A^)-va\ue 
is shown (in log-scale). Plots correspond to billiards in order as depicted in Fig.[3l 




Figure 10: (color online) Data from Fig. [9] shown for each axis separately. Left: distributions of Ap-values for 
all billiards (indicated by different symbols) done over the entire /c-interval from A; = 85 to A; = 94. Right: the 
same for A^-values. 



much less visible. 

4 Discussion and Conclusions 

We carried out a computational experiment, designed 
as a model of scattering of quantum particles through 
billiards with variable geometries, i.e. variable classical 
integrability properties. In particular, we considered a 
fully integrable regular triangular billiard, and three 



non-integrable billiards: triangular pseudo-integrable, 
triangular weak-mixing, and Sinai-type strongly chaotic 
(K-system). Billiards were designed with a double-slit 
opening on their base in a way to allow for two scat- 
tered probability currents to interfere and create an 
interference pattern observable on a screen. The prop- 
erties of the interference patterns were investigated in 
relation to the classical integrability of the billiards. 
Experiments were done for a range of particles' ener- 



11 



gies, and the interference properties of the obtained 
patterns were studied statistically, comparing among 
the billiards. 

Dramatic difference between the interference pat- 
terns exhibited by the integrable (regular) billiard and 
non-integrable billiards was found. In the integrable 
case an overwhelming fraction of the observed pat- 
terns is constructive, i.e. resembling the outcome of 
the ideal plane- wave double-slit experiment. In con- 
trast to this, non-integrable billiards display a vast 
variety of interference patterns, with a large fraction 
of them being asymmetric or of low visibility (con- 
trast ratio). We substantiated our study by introduc- 
ing the quantities Ap and in Eq. dSj) that quantify 
the phase-symmetry and intensity-symmetry between 
the two slits. As clear from Fig.[9l the integrable bil- 
liard exhibits a sharp peak in the joint distribution of 
these parameters at {Ap,A^) = (0,0), whereas non- 
integrable billiards show much more dispersed distri- 
bution of (Ap, A<^). 

Although all the experiments were conducted with 
a fixed height of the potential barrier, we investigated 
the effect of changing its height. Raising the barrier 
makes the resonant peaks (cf. Fig. [5]) narrower and 
gives a better clarity of the interference pattern, as this 
way the particle "feels" the billiard's geometry better. 
However, this also reduces the total output flux for all 
/c- values, making many interference patterns unusable 
for statistical study. Our selection of the barrier height 
optimizes between these two issues. 

The interference patterns change smoothly with k. 
The fact that non-integrable billiards display a larger 
variety of different patterns, also testifies about their 
more rapid change with increasing k. Additional in- 
sights into the mechanism of quantum dephasing might 
be obtained by studying the speed of the structural 
change of the interference patterns with k. 

Along the lines of conclusions reached in [1], we 
thus re-confirmed the connection between the integra- 
bility of a classical system, and the appearance of quan- 
tum dephasing responsible for the destruction of the 
interference patterns. The classically chaotic geome- 
try seems to be able to generate dephasing in both 
time-domain (as shown in [1]) and in energy-domain. 
Interestingly, while the distinction between integrable 
and non-integrable billiards was clearly revealed, the 
difference among the non-integrable billiards is much 
less visible. As it appears from Figs.[9]&:[10l the pres- 
ence of non-constructive interference patterns is sim- 
ilarly prominent for all three non-integrable billiards. 
It seems that all three of them are equally able to gen- 
erate dephasing, and no clear conclusion can be made 



from these results regarding the strength of the dephas- 
ing induced by billiards of different non-integrability 
levels. Further work and possibly new experiments are 
needed to clarify this interesting question. 

The computational experiment presented here can 
be easily constructed and carried out in a real labora- 
tory. 
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